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Abstract 

A variety of observations of novae are discussed in 
light of our theoretical models. We propose that the nearly 
constant boloraetric luminosity of FH Ser observed by 
Gallagher and Code originates in the non-degenerate hydrogen- 
burning region at the bottom of the hydrogen-rich envelope 
which remains after the primary ejection. The shift of the 
wavelength of peak emission from the visual to shortward of 
the ultraviolet is caused by the decrease of the photospheric 
radius of the remnant envelope as the bolometric 
luminosity stays nearly constant. The oscillations in the 
light curve of GK Per during the transition stage can be 
explained by a pulsation of the remnant envelope when it is 
the size of the Roche lobe. The CNO over-abundances in 
novae reported by various observers are strongly suggestive 

of our nova mechanism. Finally, the implications of the 

13 15 

upper limits of C and N found by Sneden and Lambert in 


DQ Her are discussed. 
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Subject Headings: Abundances - Hydrodynamics - Hass Loss - 

Novae > Nucleosynthesis - Ultraviolet. 
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I. INTRODUCTION 

* 

In a recent series of papers (Starrfield, Truran, Sparks 

and Kutter 1972; StarrfJcld, Sparks and Truran 1974a, b; 

hereafter, called papers I, II, and III, respectively) we 

have discussed our theoretical models of novae outbursts. 

Assuming the Kraft model (1963), which has a white dwarf 

accreting hydrogen-rich material from a binary companion that 

has overflowed its Roche lobe, we have followed the 

hydrodynamics of a thermonuclear runaway in the degenerate 

hydrogen-rich envelope of the white dwarf by means of an 

implicit hydrodynamic stellar evolution computer code (Kutter 

and Sparks 1972). The decay of the 8^-unstable nuclei formed 

27 29 

during the thermonuclear runaway causes 10 - 10 grams of 

the hydrogen-rich material to be ejected with velocities of 
200-2400 km s”^. These values agree with the observations 
of novae. We have also reproduced the general feafires 
in the early portion of the light curve of a fast nova 
outburst. Upon finding that a less luminous white dwarf 
requires a longer time to thermonuclear runaway and has 
a more intense outburst, we have presented a qualitative 
explanation for the Kukarkin-Parenago relationship (paper II), 
as applied to the common and recurrent novae. Lastly, we 
have explained the oval shape of the nebulosity surrounding 
DQ Her (Sparks and Starrfield 1973) by taking into account 
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the interaction between the material expanding off the white 
dwa|rf and the accretion disk surrounding it. 

In this paper we wish to extend the results from our 
hydrodynaialc models and show that they can explain some 
rather diverse observations of the nova outburst. In ( II 
we will consider the recent ultraviolet observations of 
FE Ser 1970. In § III an explanation for the oscillations 
of the light curve of GK Per during the transition stage 
will be given. The elemental and isotopic abundances of 
carbon, nitrogen and oxygen produced by a nova outburst 
will be discussed in § IV. 
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II. ULTRAVIOLET OBSERVATIONS OF FH SERPENTIS 1970 

Recent ultraviolet observations of FH Ser by Gallagher 
and Code (1974) showed that beginning 4.4 days after maximum 
the bolometric luminosity remained at approximately 

4 

1.6 X 10 Lq for 53 days of observation.*?. They also found 
that in the same time period the wavelength of peak emission 
Shifted from the visual to shortward of 2000S. Infrared 
observations by Geisel, Kleinmann and Low (1970) indicated 

4 

a peak luminosity of 5 x 10 L^ 90 days after maximum. 
Therefore, the fast decline in the visual magnitude for 
FH Ser was caused by the radiative energy shifting to other 
wavelength regions not by a decline in the total energy 
output of the nova. This may, in fact, be true for all 
novae (Gallagher and Starrfield 1976). 

Let us consider a typical light curve of one of our 
hydrodynamic nova models (No. 3 from paper II) in Figure 1 
where the latter part of the light curve has been artificially 
Smoothed to remove the coarse zoning effects which cause 
large variations in the luminosity that mask the true 
behavior (Sparks 1969). The initial rise of the bolometric 
light curve is due to a shock wave initiated by the thermo- 
nuclear runaway. The light curve rapidly drops after passage 
of the shock wave but rises again rapidly as the outer layers 

are heated by the decay of the, 6^-unstable nuclei (^^N, ^^0, 

15 17 

0, F) formed during the thermonuclear runaway and mixed 
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to the surface by convection (cf. paper I). (This rapid 
decline and rise occurs on too short a time scale to appear in 
Figure 1.) The effect of this heating is demonstrated by 
the positive luminosity gradient in the outer layers shown at 

3 

a time of 2 x 10 sec in Figure 2. (All times refer to the 

4 

time after the initial rise.) At a later time, 1.9 x 10 sec, 

all of the 0^-unstable nuclei have decayed and the outer 

escaping layers have cooled due to expansion and radiative 

losses producing the dip in the light curve also shown in 

4 

Figure 1. At 7.7 x 10 sec in Figure 2, we see that the 
high inner luminosity has progressed outward while the region 
where the zones are optically thin (to the right of the left- 
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hand bracket) has moved inward. By 2.78 x 10 sec this 

"luminosity wave" has moved through all the zones giving the 

final rise in the light curve in Figure 1 (also see Sparks 

1969). This luminosity wave is caused by non-degenerate 

hydrogen burning at the bottom of the remnant envelope . 

Much of the radiation in the initial sharp peak shown 

in Figure 1 is in the ultraviolet and the computed visual 

light curve has the Initial rise, pre-maximum halt and final 

rise structures that are observed (cf. McLaughlin 1939). 

The cause of the decline in the visual light curve can be 

seen if we plot radius versus interface number in Figure 3. 

4 5 

This plot shows that from 7.7 x 10 sec to 2.78 x 10 sec 
the photospheric radius (where t ■ 2/3 and shown by the 
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left-hand bracket) is decreasing as the ejected material 
becomes optically thin. The decrease in the photospheric 
radius and the nearly constant bolometric luminosity (produced 
by the constant nuclear burning shell source) leads to an 
Increase in the effective temperature and a decrease in the 
amount of energy appearing in the visual. In order to compare 
our results with the ultraviolet observations of Gallagher 
and Code (1974) we removed the outer e.scaping material and 
evolved only the remnant envelope. We do this since the 
velocities of the expanding nebula are so large that the 
change in the density allowed per time step forces the time 
step to be very small and thus effectively prevents further 
evolution of the model. This remnant soon evolved to the 
"equilibrium model" in Figure 3 which has a luminosity of 

4 

1.7 X 10 in agreement with their observations. 

The character of the energy production has now changed 

from fast CNO- burning in a degenerate envelope to the normal 

CNO cycle in a non-degenerate material. The burning is 

occuring in a thin shell. The energy is coming from proton 

capture rather than the decay of the S‘*’-unstable nuclei, the 

shell source temperature has declined from its peak value 

of ~2 X lO^^K to a value of -.5 x 10^°K, and the energy 

16 —1 —1 

generation rate has dropped from 10 erg gm sec to 

a 

10 erg gm sec . At this energy generation rate, it 

3 

would take -4 x 10 yr to convert all the hydrogen to helium in 


a 
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the remnant envelope (-10 if there were' no further 

mase loss. However, such a high luminosity can cause mass 
loss by means of radiation pressure (see e.g. Sparks and 
Kutter 1972; Kutter and Sparks 1974). (A proper study of 
radiation pressure mass loss requires the evolution of a finely- 
zoned model which is in progress.) In fact, because the 
bolometric luminosity of two 50 year old novae (Gallagher 
and Holm 1974) is between 10 to lOOL^ and a 7 year old nova 
is less than 500L^ (Gallagher 1974), it appears that the 
observed remnant envelope does not evolve on a nuclear time 

scale. In addition, it is important to note that it will 
be necessary for the remnant hydrogen envelope to be either 
expelled, burned or cooled before the star can return to 
conditions where another thermonuclear runaway can occur and 
produce an outburst. 

Gallagher and Code (1974) also found that the wavelength 
of peak enission shifted shortward with time. Using Weln's 
law we have derived effective temperatures for their obser- 
vations and have plotted both our theoretical and their 
observed values. in Figure 4. The rather large error bars 
for the last point of the 45-zone nova model are caused 
by the normal difference procedures which assume uniform 
temperature zones. The removal of the outer escaping material 
i-n order to evolve the. remnant envelope causes an increase in 
the initial effective temperature so that the remnant calcu- 
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Ifttion represents an upper limit on the effective temperature 
as a function of time. The final effective temperature for the 
remnant envelope (log T ■ 4.8) can also be compared with 

G 

the values for V603 Aql (log T * 4.40) and RR Pic (log T > 

6 G ' 

4.54) found by Gallagher and Holm (1974). They also point 
out that if they had compared their observations with 
spherical atmospheres, they would have found even higher 
effective temperatures. As can be seen in Figure 4, the 
effective temperature of both of the theoretical models 
Increases faster than the observed effective temperature. We 
explain this by the fact that the computer code (Kutter and 
Sparks 1972) is one-dimensional and cannot take into 
account the existence of the binary companion. Since the 
effective gravity is less for the white dwarf in a binary 
system than for a single star, it is not surprising that the 
radius of the remnant will decrease less rapidly and thus 
the effective temperature increase more slowly than in our 
models. 

Finally at 90 days after maximum visual light Geisel, 
Kleinmann, and Low (1970) observed that the infrared radiation 

due to dust at about 900^K, peaked at 5 x lO^L^. They 
suggested "the presence of a previously undetected source 
of energy comparable to or greater than the energy released 
in optical emission. . . A hydrodynamic nova model 
including grain formation is impossible at this time. 
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However, If grain formation in the nova ejecta takes place 
on a time scale of days as they suggest then it is reasonable 
to assume that the grains are reradiating the far ultraviolet 
flux from the white dwarf in the infrared. 
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III. OSCILLATIONS OF THE LIGHT CURVE OF GK PERSEl 

During the transition stage (see McLaughlin 1939) some 
novae show oscillations in their light curves. GK Per is 
the best example of this phenomenon showing quasl-periodic 
oscillations of about 4 days and an am>>litude of about 
1.5 magnitudes (Campbell 1903). These began 23 days after maxi- 
mum when the nova had declined about 4 magnitudes. It is 
the purpose of this section to demonstrate that the highly 
extended remnant envelope could be responsible for these 
oscillations when it is approximately the size of its Roche 
lobe. 

First it is necessary to calculate the pulsa^ional 
period of the remnant envelope when it is the size of the 
Roche lobe. We begin by determining the semi-major axis 
from Kepler's law 


. 43.5 (M^o ♦ Mggp) p2. 

.10 


( 1 ) 


where A ir the semi -major axis in 10 cm, P is the 
period in hours and and are the masses of the white 

dwarf and the red binary companion, respectively, in solar 
units. Two solutions for the spectroscopic orbit of 
GK Per are in the literature. Kraft (1964) found a highly 
ecceutric orbit with a period of 1.904 days while Paezynski 
(1965) found a solution with a circular orbit and a period 
of 0.685 days. Following Mumford (1967) and Gallagher and 
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Olnas (1974), we adopt Paczynskl's solution which gives 

assume that both theoretical 

(Paper III) and observational (Warner 1973) grounds 

(since it enters as the cube root its value is not critical 

to this discussion). Thus, from equation (1) we find A = 

2.6 X 10^^ cm. The equivalent Roche lobe radius for 

the white dwarf is then 1.3 x 10^^ cm (ECopal 1959, p.l36). 

The mass of the hydrogen envelope, obtained from our 

30 

•theoretical models, is approximately 2 x 10 gm which 

-4 -3 

implies a mean density of -<2 x 10 gm cm . We are not 
Including the mass of the white dwarf since it is very 
small and massive and would not participate in any pulsation. 
The pulsation period is related to the square root of the 
mean density by means of the pulsation constant Q. We use 
the value Q »• 0.052 days for Dumezil and Curien's model 
(Ledoux and Walraven 1958, p. 472) because the ratio of the 
density at the bottom of the hydrogen envelope to the mean 
density fits their density ratio the best. In fact, this 
value is close to the average so that the other Q values 
are within a factor of 2. Using this value of Q we find 
a period of 4.4 days for the remnant envelope when it is the 
size of the Roche lobe, which agrees with the observations. 

Next we compare the radius of the calculated remnant 
envelope as a function of time with the equivalent Roche 
lobe radius of GK Per in Figure 5, The remnant envelope 
radius rapidly expands past the Roche lobe radius and then 
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more gradually decreases below it. Other of our remnant 
model calculations as well as Nariai's (1974) analytical 
calculations show a similar behavior. As the remnant 
envelope expands past the Roche lobe it forms a common 
envelope about both stars since the red binary companion has 
already filled its Roche lobe. We are suggesting that the 
remnant envelope pulsates when it shrinks back to its 
Roche lobe size. The beginning of the transition stage for 
GK Per which is the time that the oscillations begin is indicated 
by the arrow in Figure 5. The theoretical calculations 
differ with the observations in two respects. First the 
time where the remnant radius crosses the Roche lobe radius 
occurs too soon compared with the beginning of the transition 
stage. Second, the time interval that the remnant radius 
is approximately the Roche lobe radius is too short. Both 
of these discrepancies are due to the neglect of the binary 
companion in the theoretical calculations. The effective 
gravity is less for a white dwarf in a binary than for a single 
star (see §11). In an actual event the remnant radius would 
decrease less rapidly and cross the Roche lobe radius later. 

As the remnant envelope shrinks past the Roche lobe it changes 
from a common envelope around both stars to an envelope 
around just the white dwarf. It would not be surprising if 
the time interval for this phenomenon were of considerable 
length. 


It is not our purpose to propose a specific mechanism 
to drive the pulsation but merely to point out that a 
pulsation of the remnant at the Roche lobe radius has 
the proper perio * It may not even be necessary to have a 
driving mechanism since this is a short lived phenomenon. 

The oscillation quasi-period of GK Per also shows a 
tendency to increase with time. This could be due to the 
decrease of the mean density as the envelope loses its mass. 
Why does GK Per show this phenomenon more clearly than any 
other nova? The answer may lie In the fact that the white 
dwarf of GK Per has a bigger Roche lobe than any other cata- 
clysmic binary except for T CrB. Unfortunately, there is no 
other nova that showed oscillations in its light curve and has 
had an orbital period determined. 
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IV. ELEMENTAL AND ISOTOPIC ABUNDANCES OF CNO 

In papers I and II we have demonstrated that an enhance- 
ment of CNO nuclei in the material undergoing a thermonuclear 
runaway is necessary for ejection. Comparing our predicted 
abundances in the ejecta with solar abundances (Cameron 
1968), we found that carbon was four times solar, nitrogen 
sixty times solar and oxygen three times solar. The isotopic 
ratios of CNO nuclei in the ejecta were also distinctive. 
Observational evidence bearing upon these features is reviewed 
in this section. 

As noted in Paper I, Pottasch (1959) observed that 
carbon was solar, nitrogen forty-five times solar, and oxygen 
five times solar for five novae, while Mustel* and Boyarchuk 
(1959) and Mustel' and Baranova (1965) found that the CNO 
nuclei were one hundred times solar for DQ Her. In 
addition, Antipova (1969) argued that a CN enhancement was 
necessary to explain the appearance of the CN absorption bands 
in DQ Her. For the recent slow nova HR Del 1967 different 
observers have found that carbon and nitrogen were ten 
times solar (Ruusalepp and Luud 1971), oxygen was greater 
than three times solar (Sanyal and Robbins 1975), nitrogen 
was five to ten times solar (Anderson and Gallagher 1975), 
and carbon was approximately ten times solar while nitrogen 
and oxygen were approximately one hundred times solar 
(Antipova 1974). Ruusalepp and Luud (1971) also studied 
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the emission spectrum of nova and planetary nebulae and 
found that oxygen and nitrogen are more abundant in novae 
than in planetaries. Finally, Ruusalepp and Luud (1971) 
found that the fast novae (and therefore the brighter novae) 
are more nitrogen enhanced than slow novae. This agrees 
with our theoretical studies, since our models with larger 
CNO enhancements were brighter, faster, and necessarily more 
energetic. The range in values determined for HR Del measures 
the difficulties associated with such observational efforts. 
Nevertheless, we feel that the large number of observations 
of CNO over-abundances in novae provide strong support for 
our theoretical models . 

Another of our theoretical predictions is that the con- 
X3 15 17 

centrations of C, N, and 0 will be greatly enhanced in 

nova ejecta. An average over our models reveals: 74% of 
13 15 

the carbon is C, 61% of the nitrogen is N, and 21% of the 
17 

oxygen is 0. These correspond to isotopic enhancements 
relative to solar material by factors of 66, 160, and 570, 
respectively. Such isotopic anomalies are, in fact, an 
tinambigious signature of CNO hydrogen burning at high temper- 
atures on hydrodynamic time scales, and provide a potentially 
powerful test of our theoretical models. Unfortunately, there 
are major difficulties associated with obtaining isotopic 
abundances in nova ejecta. Sneden and Lambert (1975) have 
recently found that ^^C/C < 40% or ^®N/N < 33% or ^^C/C <: 25% 
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and N/N < 25% for the very slow nova DQ Her. They obtained 
these values from analyses of the CN molecular band at 3883^ 
which appeared in the spectrum near maximum light. 

In their concluding remarks Sneden and Lambert interpret 

their upper limits as approximate estimates of the isotopic 

abundances and on this basis find them in disagreement with our 

models. The point of disagreement arises from an investigation 

of the effects of accretion on the progress of the outburst 

(Starrfield, Sparks and Truran 1974c). In this study, we 

Included the effects of only the recently accreted material 

in the calculations and neglected all of the fine details of 

the accretion process (e.g., the emitted spectrum and the 

accretion shock structure). As a result of the infall, a 

temperature inversion was produced in the white dwarf envelope 

which effectively prevented the CNO enhanced material from 

reaching the surface layers. Since there is no mixing during 

the later stages of the expansion, we found that the surface 

layers showed no evidence of our ejection mechanism - enhanced 
13 ^5 17 

C , N and 0 . This is in apparent disagreement with the 

findings of Sneden and Lambert if we treat their numbers, 
as they did, as firm values. 

In an actual event, the temperature inversion would 
probably not be realized when accretion occurs over a long 
time period and CNO enhanced material would probably reach 
the surface by maximum light. It might appear that there 
is now a discrepancy in the opposite sense — our calculations 
predict too much of an isotopic enhancement. We believe that 


this discrepancy is not very serious for the following 
reasons. Very accurate descriptions of the accretion process 
and the mixing process (i.e., time-dependent convection and 
overshooting) will be necessary to determine the exact 
amount of processed CNO material seen at maximum light. In 
addition, the CNO enhanced material coming off the white 
dwarf will collide and mix with material in the ring surrounding 
the white dwarf, thus reducing the enhancement. This possibility 
is especially applicable for DQ Her since the inclination of 
the binary orbit (and thus the ring) is near 90° (Walker 1956; 
Mustel* and Boyarchuk 1970). In fact, the collision with 
the ring in this particular orientation has been utilized to 
explain the oval shape of the nebulosity now surrounding 
DQ Her (Sparks and Starrfield 1973). Finally, DQ Her was 
a very slow nova and our models seem most applicable to the 
very fast novae such as CP Pup, V 603 Agl , or Nova Cygni 
1975. We should keep in mind that while the observed upper 
limits on the isotopic abundances are a factor of 2 or 3 

lower than our theoretical values (based on models without 

13 15 

accretion), C/C is a factor of 36 larger and N/N is a 

factor of 89 larger than solar system values. Also, normal 
CNO burning will produce only ^^C/C ^.20 and ^®N/N “ 2-4 x 
(Caughlan and Fowler (1962). 

It is clear from the abundance results presented in this 
section that novae show strong evidence of the presence of the 
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enhanced concentrations of CNO nuclei that we have found 
necessary to produce a nova outburst. Unfortunately, 
this does not constitute proof that these nuclei have been 
processed through a high temperature CNO hydrogen burning 
region. Such proof awaits a definitive study of CKO 
isotopic abundances in the nova ejecta. 
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V. SUMMARY 

We have used the results of our hydrodynamic nova models 
to explain a number of diverse observations. The nearly 
constant bolometric luminosity of FH Ser observed by 
Gallagher and Code (1974) originates within the non-degenerate 
hydrogen-burning shell source at the bottom of the hydrogen- 
rich envelope which remains after the primary ejection. The 
shift of the wavelength of peak emission to the ultraviolet 
is due to the decrease of the photospheric radius of the 
remnant envelope as the bolometric luminosity stays nearly 
constant. The oscillations in the light curve of GK Per 
during the transition stage can be explained by a pulsation 
of the remnant envelope when it is the size of the Roche 
lobe. The CNO over-abundances in novae reported by various 
observers are strongly suggestive that our nova mechanism 
is actually operating in real novae. Nevertheless, we shall 
have to wait for a definitive study of the isotopic abundances 
before we can be certain that high temperature CNO hydrogen 
burning has occurred during the outburst. 
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Fig. 5 


FIGURE CAPTIONS 

. — The computed bolometric and visual light curves of 
a hydrodynamic nova model (No. 3 from paper II). 

.^Luminosity as a function of interface number and 
time. The time refers to the time since the 
initial rise and the left-hand bracket indicates 
where the optical depth is 2/3. 

. — Radius as a function of interface number and time. 

The time refers to the time since the initial rise 
and the left-hand bracket indicates where the optical 
depth is 2/3. The "equilibrium model" is the remnant 
after it has reached nearly hydrostatic and thermal 
equilibrium. 

. — Effective temperature versus time. The squares 
represent observed values from Gallagher and Code 
(1974) assumming Wein's law. The triangles and 
the circles represent the calculated values from 
the nova model and the remnant model, respectively. 

.— The outer radius of the calculated remnant model 
as a function of time. The dashed line represents 
the equivalent Roche lobe radius of Gi; Per and the 
arrow represents the beginning of its transition 


stage. 
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